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Background 


• Atomic force microscopy (AFM) has rapidly evolved from 
a quasi-static method for assessing material properties 
at the nanoscale into a fully dynamic technology 

• Various dynamic implementations of AFM, collectively 
called dynamic AFM or d-AFM, have been developed to 
assess surface, near surface, and subsurface properties 
of the material such as elastic moduli, adhesion, 
viscoelasticity, embedded particle distributions, device 
integrity, and topography 
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Equipment arrangement for various d-AFM 
techniques using either one or two probing signals 
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Cantilever-sample interaction force 



• The d-AFM output signal is derived from the interaction 
force between the cantilever tip and the sample surface 
that is dependent on the tip-sample separation distance z. 


• A typical force-separation curve (obtained from Lennard- 
Jones potential) is shown below 
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d-AFM image degradation 

Certain conditions of d-AFM operation generate chaotic 
oscillations that result in image streaking, as shown below 
for a RDF-AFUM micrograph of 50mm-thick single wall 
carbon nanotubes embedded in a LaRC-CP2 polyimide 


matrix. 
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Bi-stability 



The existence of two stable states of oscillation at certain 
frequencies (bi-stability) has been proposed as a possible 
cause of chaotic oscillations, resulting from interaction force 


nonlinearity. 
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Research objective 


• Develop broad-ranging explanation for chaotic 
oscillations and image degradation in d-AFM 
micrographs based on solutions to cantilever 
dynamical equation that includes quadratic 
nonlinearity in the interaction force. 

• Derive dynamical conditions under which 
subharmonic generation leading to chaotic 
oscillations is assured. 
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Cantilever dynamical equations 

The cantilever dynamical motion may be expressed quite 
generally as an expansion over vibrational modes n of the 
form 

y(x,t) = 2 Jn= Y n (x)rj n (t) (i) 

where 

x = spatial coordinate t = time 

Y n (x) = orthogonal basis set of cantilever shape-dependent 
spatial eigenfunctions 

r] n (t) = temporal components of the vibrational modes 
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Cantilever dynamical equations 

The temporal components rj n of the vibrational modes are 
obtained from the solutions of 



in 


- + Y„ — — + kj i = F cos (Ot + F 

dt 2 dt 


dz\ 
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( 2 ) 


where 

m c = effective mass of cantilever y n = damping coefficient 

k n = cantilever free-space stiffness constant 

F c = cantilever driving force z = cantilever-sample 

separation 

F(z,dz/dt) = cantilever-sample interaction force 

Note: dz/dt accounts for the transfer of energy between the 
cantilever and the sample surface. 
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Interaction force expansion 

Let z 0 = the set-point (rest) position of the cantilever tip on 
the force curve and write z- z 0 = rj n for mode n. 



Expand the interaction force F(z,dz/dt) as 
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FjQ (/ = 1, 2, 3, ■■■) = ‘spring stiffness’ coefficients 
F 0 j (j = 1 , 2, 3, •■■) = coefficients that account for the transfer 
of energy from the cantilever to the sample. 
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Interaction force expansion 

The interaction force expansion is not a Taylor series. 

Taylor series is limited to small cantilever displacements. 

More efficient expansion is obtained by placing interaction 
force function in the space of smooth functions over entire 
range of cantilever displacements, equip the space with an 
inner product defined for arbitrary functions g and h as 

(g,h) = fr r g(ri,fi)h(ri,fj)dridfj ( 4 ) 

^7 max ^7 max 

where rj max and r7 max dot are the maximum cantilever 
displacement and the maximum rate of change of 
displacement, respectively. 
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Expansion coefficients 

To obtain Fy, we form the inner product <(F-f),(F-f)> where 
F(ri,drj/dt) is the actual function and f is the polynomial 
expansion of F truncated to second order, take derivatives 
with respect to each of the coefficients Fy in the expansion, 
set each derivative to zero, and solve the resulting system 
of linear equations. 



For present purposes the most important coefficient is F 20 
given as 


F = 

1 20 


15 


[3(F,ri 2 )-vL(F,l) 



Previously shown that F 20 leads directly to bi-stability. 
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Dynamical equation 

Substitute polynomial expansion to second order in the 
temporal equation for vibrational mode n, consider only 
fundamental mode n = 1 , drop the subscript, and divide by 
m c to obtain the dynamical equation 

+ T— + (o^rj - /3r/ 2 = F m cos cot (6) 

dt dt 

where 

r = {Y-F 0 i) /m c <o\ = {k - F 10 ) / m c F m = F c /tn c 

fi = F 20 / m c 
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The Melnikov method 


• Under certain conditions, solutions to the dynamical equation include 
subharmonics and chaotic oscillations. 


• The Melnikov method provides means to determine threshold 
conditions necessary to assure a cascade of period doubling 
bifurcations (subharmonics) leading to chaos in a perturbed nonlinear 
system by assessing the distance between stable and unstable paths 
in phase space (velocity-position plots) emanating from an unstable 
limit cycle (if the limit cycle exists). 



position 
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Threshold conditions for chaos 


• Detailed application of the Melnikov method to the 
dynamical equation to be published elsewhere. 

• Most important result is a quantitative assessment of the 
threshold drive amplitude that assures a cascade of 
frequency-halving bifurcations (i.e., subharmonics coll, 
co/4, co/8, of the cantilever drive frequency co) 
terminating in chaos. 


• The threshold cantilever drive force per unit mass ( F m ) th 
given as 
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Note dependence on parameter of nonlinearity /3. 
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rj max > Wdt) max dependence of (FJ th 

• Threshold condition {F m ) th = (r<yf / 5jto> 2 P)$,mh(<x>ji / (o x ) assures the 
generation of a cascade of period doubling bifurcations leading to 
chaos as drive amplitude F m monotonically increased above ( F m ) th . 



• (F m ) th depends on the cantilever drive frequency <x>, the natural 
oscillation frequency co v the damping/energy transfer coefficient r, 
and the nonlinearity parameter of the interaction force. 

• Recall P = F 20 /m c and ^>0 = ( 15/16j 7max J 7max)[ 3 ( F > J 7 2 )- J 7max( F > 1 )] ■ 

5 

• Conclude (F m ) th depends on maximum displacement of 
the cantilever from its set-point and the maximum rate of 
change of the displacement over oscillatory cycle. 
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Conditions for scan instability 

Quantitative assessment of threshold drive force (F c ) th = m c (F r 
be calculated from the Lennard-Jones potential 
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where 0 O = 4.05 x 1 0 22 J and z m = 0.27 nm. 

The interaction force F(z,dz/dt) = F^d^/dt) = - (d0 LJ /dz). 

Substitute F(^,d^/dt) in the expression for F 20 
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to obtain j5 = F 20 lm c 

Calculate ( F c ) th from threshold equation 


( F c),h =m c{ F m) th =m c {T(oll5n(o 2 p)smh{(onl(o l ) 
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Log-linear graph of\(FJ th \ versus setpoint z 0 
for oscillation displacement amplitude 

hmax ~ 0.1 nm. 
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Log-linear graph of (FJ th versus oscillation 
displacement amplitude rj max for fixed set-point 

z 0 = 1.0 nm 
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RDF-AFUM image degradation 

RDF-AFUM micrograph of 50mm-thick single wall carbon 
nanotubes embedded in a LaRC-CP2 polyimide matrix. 
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Chaos in d-AFM 

• Present model applies to all commonly used d - AFM modalities. 

• For example, contact models of cantilever dynamics show that 
chaotic oscillations occur when the cantilever makes contact with a 
hard sample surface - a problem of special importance for AM-AFM 
(tapping mode). 

• Present model based on shape of the force curve as a continuous 
curve at the atomic scale, 

• ‘Contact’ in model occurs when the cantilever set-point is in the 
more linear region of the interaction force curve but the amplitude of 
cantilever oscillation is sufficiently large to produce an excursion into 
the highly nonlinear region in the vicinity of the ample surface and 
the minimum of the curve. 

• The large amplitude excursion produces a dramatic increase in the 
value of /5 and a decrease in (F m ) th that results in chaotic oscillations 
in AM-AFM. 
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Conclusion 

• Developed explanation for chaotic oscillations and image 
degradation in d-AFM micrographs based on solutions to 
cantilever dynamical equation that includes quadratic 
nonlinearity in the cantilever-sample interaction force. 

• Represented interaction force as polynomial expansion 
with coefficients that account for the interaction stiffness 
parameter, the cantilever-to-sample energy transfer, and 
the amplitude of cantilever oscillation. 

• Derived dynamical conditions under which period 
doubling bifurcations leading to chaos in d-AFM is 
assured. 
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